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Abstract 

We consider the partition function of the superconformal Chern-Simons theories with 
the quiver diagram being the affine D-type Dynkin diagram. Rewriting the partition 
function into that of a Fermi gas system, we show that the perturbative expansions in 
1/N are summed up to an Airy function, as in the ABJM theory or more generally the 
theories of the affine A-type quiver. As a corollary, this provides a proof for the previous 
proposal in the large N limit. For special values of the Chern-Simons levels, we further 
identify three species of the membrane instantons and also conjecture an exact expression 
of the overall constant, which corresponds to the constant map in the topological string 
theory. 
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1 Introduction 

The ADE classification not only provides a beautiful structure in the mathematical science, 
but also plays important roles in M-theory. As an example related to the M5-brane, there is 
an ADE classihcation of six dimensional M = (2, 0) theories, coming from the classihcation 
of the singularities on which the M5-branes are placed [1]. Also, even for the six dimensional 
M = (1,0) theories, the ADE classification continues to be crucial [2]. 

In the context of the M2-branes, the three dimensional U{N) quiver gauge theories can also 
be classified by the affine ADE-type, or ADE-type, Dynkin diagrams. A large class of three 
dimensional superconformal Chern-Simons theories can be constructed by quiver diagrams as 
follows. For each vertex in the quiver diagram, we assign a vector multiplet of a gauge group 
U{N), with N proportional to the comark. For each edge connecting two vertices, we assign a 
pair of hypermultiplets which are in the bifundamental representation under the gauge groups 
on the two vertices. For example, the simplest quiver, Ai, gives the ABJM theory [3]. For 
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these theories, the localization technique allows us to express the partition function on as a 
hnite dimensional matrix model [1]. Then, it was shown in [5] that, if we require that the long 
range force among the eigenvalues vanishes (also known as the balance condition in [6]), the 
quivers have to be of ADE-tjpe. The behaviour of these matrix models in the limit N ^ oo 
was studied for the A-type quivers in mi, for the D-type quivers in [Sl[ni[in] and for the 
T'-type quivers in [10]. Interestingly, they observed a universal scaling law in the partition 
function for all the ADE-type quivers 


Z{N) ~ exp 



( 1 . 1 ) 


with C some constant depending on the quivers. This scaling law is a characteristic property 
of the M2-branes in the context of the AdS/CFT correspondence [TT] . 

For the ABJM theory, after this leading N 2 scaling behaviour was obtained in [12], the 

3 

partition function was studied in full detail. It was shown that the N 2 behaviour is completed 
by the perturbative 1/A^ corrections into an Airy function |13] . 


= e^C-^ Ai\c-^{N - B) 


( 1 . 2 ) 


with some constants A [IT], B and C, up to non-perturbative corrections in 1/N. Later it 
was pointed out that the partition function of this theory can be rewritten as that of an 
ideal Fermi gas system with N particles [15]. This formalism is so efficient that it not only 
provided a simple rederivation of the above result, but even allowed the exact analysis of the 
non-perturbative corrections [T6H2T]. 


In [15] the authors also showed that the Fermi gas formalism works for theories of general 
A-type quivers, and that the completion by an Airy function is universal for these theories. 
The superconformal Chern-Simons theories of the A-type quivers were also studied in detail, 
uding the perturbative coefficients A, 5, C* and various non-perturbative corrections [sa¬ 


me 


0 


It was further conjectured in [TS] that the completion by an Airy function is universal even 
for other theories of the M2-branes. However, in the case of the DE-type quivers, it was not 
trivial whether the universal behaviour of the Airy function is valid because of the lack of the 
Fermi gas formalism due to its non-circular structure and the non-uniform comarks. 

In this paper we consider general H-type quivers, and obtain a positive answer to this 
question. Here we shall explain our setup. We consider the quiver superconformal Chern- 
Simons theory whose quiver is the Dr Dynkin diagram. We call the two vertices on the left 


* The Nf flavor matrix model also allows the Fermi gas formalism which turns out to be equivalent to that 
for some of the A-type quiver. These matrix models were studied in [261128] . The equivalence between these 
matrix models is proven systematically in [29]. 
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end as jj, and /i', those on the right end as u and v', and label the others as 1, 2, • • • , r — 3, as 
in hgurelU In addition to the nnmber Nm of the gange gronp U{Nm), the Chern-Sinions level 
is assigned on each vertex, which we take as (m = 1, 2, • • • , r — 3) 

= fc(Si S 2 ), km — k(^Sm+l Sm+2}i ky = /c(Sr-_i Sr)-, 

k'^ = k{-Si - S2), kl = k{Sr-l + Sr), (1.3) 

with Sm extended to arbitrary real numbers. This choice is the most general one under the 
requirement of the superconformal symmetry of the theory j3]. Below, we shall hrst present a 



Figure 1: The label of the vertices of the Dr quiver. The number on each vertex is the rank 
of the gauge group, which is proportional to the comark. 


Fermi gas description for the partition function of this theory in section [2l In the Fermi gas 
formalism, it is rather convenient to introduce the chemical potential p dual to N and study 
the grand potential J(/i) dehned by 


E e^^Z(iV). 

N=0 


(1.4) 


Other than Sm, the grand potential is controlled by two parameters, the chemical potential /i 
and the overall Chern-Simons level k which plays the role of the Planck constant of this quan¬ 
tum statistical system, h = 2nk. The grand potential turns out to be a cubic polynomial of fi 
if we neglect the non-perturbative effects. In this manner, all order perturbative corrections 
to the partition function in 1/N are taken into account. As a result, we obtain the expression 
of an Airy function (11.21) in section [3l The coefficient C relevant to the leading N 2 behaviour 
is obtained as 


C 




r—1 

E 

m=l 


'^m+1 


GrCr+l / 


(1.5) 


where the variables a are given by 

r r 

^ ^ (I '^n| “l~ \^m '^n|) ^0 2(r 2), <^r+l 2 ^ ^ (^*6) 

n=l n=l 
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with the reordered s- 


m? 


0 < |Sr| < |Sr-l| < • • • < |Sl 


(1.7) 


This coincides perfectly with the previous proposal in ra, where the authors further tried 
to give a Fermi surface interpretation to their proposal. In the limit of /c —)■ 0, we can also 
compute the coefficient B, and the result is 


B 




+ 0{k). 


( 1 . 8 ) 


If we further restrict ourselves to the special values of Sm 

Si = S2 = ■ ■ ■ = Sr = I, (1.9) 


we can compute the higher order corrections in k to the grand potential, including the per¬ 
turbative coefhcients and the non-perturbative corrections. We analyse this model in section 
[Hand conclude that, for the non-perturbative corrections, there are three kinds of membrane 
instantons with exponents e and e“ 2 . For the coefficient A, we observe that the 

result is consistent with 


A 


^ABJM(2r/c) -|- r^AABJM(2(r — 2)k) , 


( 1 . 10 ) 


up to 0{k^), where Aabjm(^) denotes that coefhcient for the ABJM theory. Both the co¬ 
efficient A and the membrane instantons are reminiscent of the expressions for the A/" = 4 
supersymmetric theories of the A-type quiver [2^124] . 


2 Fermi gas formalism 


In this section, we shall present a Fermi gas description for the superconformal Chern-Simons 
theories of the H-type quiver. 

The partition function of this theory is given by 


Z{N) = 


^ F)27Vy(m) y 

Ivi m m ivT (2Ar)! h’ 


with the integration measure 


D/Ui 


27T 



271 


ikm 

g i-K / 


dvi ikjy 

Uui = —e 4^ 
27r 




( 2 . 1 ) 
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2ti 


Du' 


du' K 

- -p 47r 

27r 




( 2 . 2 ) 


Here the numerator V, coming from the vector multiplets in the adjoint representation, is 
given by 


N 


V = ]^f2 sinh 


fli fXj 


N 


]^f2sinh 


Hi - Hj 


r-3 2N 


n n (2 sinh 


i<j 


X 


i<j 


m=l a<b 


. 2 


N 


(2 sinh 


Ui -Uj\^ 


N 


( 2 sinh 


u' - z/' \ 2 


(2.3) 


i<j 


i<j 


and the denominator H, coming from the hypermultiplets in the bifundamental representation, 

is 


N,2N (1) N,2N 

Hi ^o, '' 


hf = ^2 cosh 


(2 cosh 


/ \(1) r-A 2N,2N 

Hi — 


2,a 


2,a 


^2 cosh 

m=l a,h 


^(m) _ ^(m+1) 


X 


2N,N 

n( 


2 cosh 


xh-3) _ 2V,iV 

'Aa ^-3 


(2 cosh 


\h-3) / 

Aa — Ui 


(2.4) 


2.1 Density matrix from matrix model 


To express the partition function fl2.1l) of the superconformal Chern-Simons matrix model of 
the Zi)-type quiver in terms of that of a Fermi gas system, in this subsection let us hrst rewrite 
the generating function of the matrix model into a Fredholm determinant. 


First, we rewrite the integrand of the matrix 


model fl2.ll) intcUl 



Here we have introduced the combined variables (/l)a=i,..., 27 v = {{H)a=i,---,Ni{H')a-N=i,-,N) 
and {u)a=i,-, 2 N = {{y)a=i,-,N){T^')a-N=i,-,N)- Namely, the second factor in (12. 5 p is the de¬ 
terminant of a matrix which is a vertical array of two N x 2N rectangular matrices, one with 

tWe knew of a related work m from the reference list of [29]. In a seminar by Nadav Drukker at Nagoya 
university, we learned that actually they had a similar idea in rewriting the integration measure into a deter¬ 
minant of hyperbolic cosecant functions as in (12.511 . We are grateful to Nadav Drukker for valuable discussions. 

^ Since the singularities appearing in the first and last determinants are originally absent in (12.4p , we expect 
them to be harmless. 
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/ M'—1 / 

components (2 cosh 2 *’ ) other with components (2 cosh " 2 ^ ) ■ Similarly, the 

second last factor is the determinant of a matrix which is a horizontal array of two 2N x N 
rectangular matrices, ((2cosh ——) ((2cosh-2-^—Note that in 

the above rewriting we have used the formulae 


>2NxN' 


n.<,2smhi^n,<.2smh=^ 


n 2 coshi!i^ 


= det 


nf<,2sinh^nf<,2sinh^ 


2 cosh 


fMi-Uj ’ 


n„-2sinh^^ 


= (_l)^(^-i)^det 


Lt,J - “***** 2 

which follow from the Cauchy identity 

- Xj) U^<Jiyi - Vj) 


2sinh 


fii-Uj > 


(2.6) 


= det 


(2.7) 


Then, using the formula proved in appendix A of [32] with the same ranks, we can combine 
the series of 2N x 2N determinants into 


U^,jixi + yj) Xi + Vj- 

by the substitutions x, = and i/j = or Uj = — 


Z{N) = 




m m m m 

Here the functions M, N and L denot 


det p'.) det vi) det u'A 


( 2 . 8 ) 




2 sinh 

r—3 


2 m=l 


2 cosh ^ 


2 sinh 


r r—4 

n 




2 cosh 




2 cosh 


a0;;-3)^’ 


(2.9) 


with the matrix in the second determinant in 02.81) given explicitly by 

Lijla^Vb) = 




( 2 . 10 ) 


Furthermore, if we use the formula in [32] for different ranks, we can perform the p' and u' 
integrations to hnd 


Z(Af) = (-l) 


N 


IV! W 


det 


L{fii,Uj) {L» N){fii,Uj) 

(M • L)(/ij, izj) {M • L» 


( 2 . 11 ) 


§The rank N of the gauge group should not be confused with the function N{v,v'). Also, the chemical 
potential ^ appearing later should not be confused with the integration variables Hi and /i'. 
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where • stands for either the /i' integration or the v' integration in fl2.2p . It should be clear 
from the context which integration it stands for. For example, (M • L • N){^i,i/j) in the 
lower-right block denotes 

{M»L» iV)(/ii, Uj) = j Dn'Du'M{ni,n')L{n', uj). (2.12) 

Now we can apply the formula flA.ip in our appendix [A] and hnd 

= p\ ( 2 - 13 ) 

with the four N x N blocks given by 

Fii = —(L oN»L + L»No fij), P12 = {LoN»L + L»NoL)» M{^i, fij), 

P21 = —M •{LoN»LpL»No L){jjii^ fij), P22 = M»{LoN»L + L»NoL)» M{jj,i, Hj). 

(2.14) 

Here o denotes the u integration. If we further introduce the chemical potential /i and dehne 
the grand potential J(/i) as (11.41) . we hnd 

= \j det(/ + e^p), (2.15) 

by using the formula in appendix [Bl Here the density matrix p is 

p = nP, (2.16) 

and the other matrices are given in flB.2p . Then it is not difficult to observe that the density 
matrix can be put into 

- _ f f f Nyy! 

Here we have regarded the functions M, N and L as matrices, and contracted the adjacent 
indices by integrations fl2.2p . without displaying • or o explicitly. After suitable similarity 
transformations and rearrangements, we can further put the density matrix into 

__ ( Ly^i\ ( 0 M^^>\ ( L \ ^ 0 

Ly/^/J 0 J \^L^Iy L y' J \jSfy'y 


Nyy>\ 

0 y 
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2.2 Operator formalism for density matrix 


In the previous subsection, we have reduced the study of the superconformal Chern-Sinions 
matrix model of the D-type quiver into an integration kernel p. To express the superconformal 
Chern-Simons matrix model in terms of a Fermi gas system with N particles, we need to 
further rewrite the integration kernel in the operator formalism. For this purpose, we introduce 
the canonical coordinate/momentum operators g and p obeying the canonical commutation 
relation 


Kfl = iH, 


( 2 . 19 ) 


with the Planck constant given hy h = 27rk and normalize the coordinate eigenstate as 


{q\q') = 2Tr6{q-q'). 


( 2 . 20 ) 


It is not difficult to spell out each block in J) fl2.18p explicitly in terms of the hyperbolic 
functions and the integrations. After rescaling the integration variables by 1/k, we hnd that 
the density matrix is given by 


P = 


(z/|p+|l/) (z/|p_|z/') 


( 2 . 21 ) 


with the operators p± 
P+ = 

2 cosh I 




2 cosh I 


-(.-iK.i^r-2-Sr-3)q^ . . . g-^p3-S2)q^ 


2 cosh I 


! tanh - 


X I _2 g-^(si+32)q= g-^(s2+si)^ 


tanh - 


,2 


X 


1 _ ^pM^2-S3)q^ . . . p^(Sr-2-3r-l)92_ ^^ pM^r-l+Sr)^ 2 ^{Sr-Sr-l)q^ 


, tanh - 


2 cosh I 


2 cosh I 


P- = 


1 _ p-4ni^r-l-Sr-2)iP _1 


-±{Sr-2-Sr-3)q^ . . . ^-^(^3-52)5^ 


2 

■2 1 


2 cosh 1 

2 cosh 1 


-j- p Tl Vt P 

g-^Pl+32)q^ , g-^P2+«l)g^ 2 

V 2 

2 

X ^ pJK(^2-S3)q^ . 

. . piK(^r-2-Sr-l)q^ t ^(^Sr-l-Sr)q^ 


2 cosh I 


tanh - 


2 cosh I 


2 cosh I 




. ( 2 . 22 ) 


In the derivation we have used the formulae 


/ I 1 I /\ 1 1 / I tanh I i 1 

(^1:^- r^\^) = T- -(g| o \^) = T 


2 cosh I 


k 2 cosh 


k 2 sinh 


(2.23) 
























Using this result we can simplify the grand potential fl2.15p . In fl2.15p the determinant 
is taken simultaneously over the functional space (or, in the operator formalism, the phase 
space) and over the two dimensional space. However, since the left two components and the 
right two components in fl2.2ip are identical, the determinant over the two dimensional space 
can be taken trivially 

+ e^^p), (2.24) 

with the density matrix purely in the phase space given by p = +p_. After performing the 

similarity transformation to move on the right end of p± in (I2.22p into the left end, 

we can use the relation 

= F(p - s^, (2.25) 


to rewrite the density matrix into 


P 


1 


2 cosh ^ 
1 

^ 2 cosh 


1 tanh ^ + tanh ^ 

2 cosh ™ 2 

1 tanh ™ + tanh M 

2 cosh 2 


(2.26) 


3 Large N behaviour from Fermi surface analysis 

In the previous section, after switching from the partition function Z{N) to the grand potential 
J(/i), we hnd that J(p) is expressed in terms of the Fredholm determinant fl2.24p of the density 
matrix fl2.26p . Since the relation is very similar to the case of the A-type quiver [T5], we expect 
that the perturbative corrections to the partition function again sum up to an Airy function 
as fll.2p with some constants C, B and A. Indeed, the expression of the Airy function follows 
from the large E behaviour of the number of states with energy smaller than E 

(3.1) 

In this section we shall show this relation, with explicit expressions of C and B up to 0{k), 
by the technique used in [T5] . 

For this purpose, let us consider the classical limit of the Fermi gas system {h 0). 
Classically, the number of states n{E) is given by the phase space volume 

niE) = ^vol{(g,p) e M^| logpo^ < U}, (3.2) 

^ The overall factor 2 compared with the case of the A-type quivers is due to the square-root in (I2.24|) . 


n{E) =tT[9{E= 2(cE‘^ + 
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where po is the classical density matrix obtained from p fl2.26p by neglecting the commntators 


log(po(g,p)) ^ = 


m=l 


log(2 cosh 


P- s^q 


+ log (2 cosh 


P + SmQ 


log(2 sinhp)^. (3.3) 


Since the right-hand side is independent of the signs and the ordering of Sm, in this section 
we assnme > 0 and fll.71) withont loss of generality. The classical Fermi snrface (i.e. the 
bonndary of the region with logpg^ < E) is plotted in hgnredl Since the region is symmetric 


P 



Fignre 2: The classical Fermi snrface for = {1,2,3} at E = 10 (the solid red line) and 
the polygon to which the Fermi snrface approaches in the limit of F' ^ oo (the dashed bine 
line). 


nnder the reflections q i—)■ —q and p i—)■ —p, below we consider only the snbregion in the hrst 
qnadrant M>q. We can fnrther divide the volnme vol{(g,p) G R>Q|logpQ^ < E} into the 
leading contribntion in the limit of large E, Vpoi, and the deviation from it, 6V, as 

n{E) = ^{Vpoi-SV). (3.4) 

Tin 

The main contribntion Fpoi can be compnted by approximating the hyperbolic fnnctions by 
rational fnnctions 


= vol< (g,p) G 


^>0 


E 

m=l *- 


\V ^mQ\ ^ |P “1“ 


2|p| < E 


(3.5) 


Since the above snbregion on the hrst qnadrant is a polygon (see hgnre[2]), whose vertices are 
located at 


2ESm \ 

’ <70 / ’ V<7m’ (Jra / 



(3.6) 
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with a given by as in fll.bp . the volume of this subregion is 

+ y (3.7) 

V^O^l rn—l ^r^r+1 / 

Now we consider the deviation from the volume of the limit polygon, 6V. First we divide 
the region between the classical Fermi surface and the polygon into the pieces around each 
line of p = Sm<i and p = 0, as in figure [3], and call the volume of each piece Vm and n^+i 
respectively, 



Figure 3: The region contributing to Vm is the region surrounded by the classical Fermi 
surface (the solid red line), the boundary of the polygon (the dashed blue line) and the solid 
black lines next to the line of p = SmQ- Each solid black line connects the origin and the 
midpoint on the edge of the polygon. 


r 

5V ='^Vm + Vr+1. 

m=l 


(3.8) 


Assuming that q and p are of order E on the classical Fermi surface, we can approximate its 
segment near the line p = s^q as 




\P 


siq\ ^ 


siq\ 


+ log 2 cosh 


P-Smq , \P + Smq\ 


+ 


- 2|p| = E, 


(3.9) 


if we neglect the non-perturbative 0{e corrections. To calculate Vm, it is convenient to 
introduce the tilted coordinate (g, p) = {q,p — Smq), 


rp+ 


Vm. — 


Mqip) - q\^), 


(3.10) 
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where q'(p) and q'{^ are the ^-coordinate of a point [q,^ on the limiting polygon of the Fermi 
surface and that on the approximant fl3.9l) . and p± are the midpoints on the edges of the 
polygon. Noting p — s^q < 0 for £ < m and p — Siq > 0 for £ > m in this piece, we can compute 
Vm as 

Vm= [ dp-^(log2cosh|-~ (3.11) 

Here, although originally the integral interval, [p_, p+] is hnite, we can replace it with the whole 
real axis (— 00 , 00 ) without affecting the perturbative behaviour in fl3.ip . since the integrand 
is exponentially small for large p. The contribution from the piece around p = 0 can be 
calculated similarly, 

TT^ 

Vr+l = - -. (3.12) 

oar+l 

Substituting these results fl3.7p . (13.Sp . (13.lip . (I3.12p into (13.4p . we obtain the large E expression 
of n{E) (13.ip with C and B given by (II.5p and (II.8p . 

So far we have been neglecting the quantum corrections. Though it is difficult to take 
care of them due to the variety of arguments of the hyperbolic functions in the density matrix 
(I2.26p . we can make the following estimation. There are two kinds of h-corrections, the 
Wigner transformation of each operator and the commutators of operators coming from the 
Baker-Campbell-Hausdorff formula. According to the Wigner transformation formula, the 
former corrections always start with the second derivatives of each term in (13.3p . Also, if the 
Hamiltonian is hermitian, since there are only nested commutators, the latter corrections again 
start with the second derivatives§ Therefore, since the second derivatives of the hyperbolic 
functions are always exponentially suppressed, the quantum corrections never change the 
asymptotic polygon of the Fermi surface in the limit of A' —?■ 00 . This ensures the behaviour 
of n{E) (13. Ij) with C uncorrected, and therefore that the perturbative partition function is 
given as an Airy function even with all order quantum corrections. On the other hand, B is 
possibly corrected due to the quantum effect. 


4 A and instantons for special quivers 

In this section, we restrict ourselves to the cases where the Chern-Simons levels are given 
by (11.31) with a uniform value of Sm (11.91) . where we set that value to be 1, which is always 
possible by the redehnition of k. Under this restriction the exact large p expansion of the 

II The requirement of hermiticity is essential also in the discussion in the A-type quiver [151123] . 
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grand potential can be computed systematically order by order in k. As a result, we obtain 
the constant part A (which appear in the partition function as fl 1.21) 1 and the non-perturbative 
corrections {0{e~^)) in the grand potential. These are just what we did in the theory of the 
A-type quiver with J\f = 4 supersymmetry enhancement [T^[23ll24] . 

For the simplicity of explanation, let us define the Hamiltonian, e~^ = 'p. After a similarity 
transformation, the Hamiltonian is given explicitly by 




(4.1) 


where we have introduced new variables 


U = log 2 cosh 


P 

T = log 2 cosh—, 


^ tanh § + tanh ^ 

S = \og - 


(4.2) 


with Q = p + q and P = p — q. Note that the Planck constant is doubled in the new canonical 
variables, [Q,P] = i{2h). 


Although we are interested in the large p expansion of the grand potential, the original 
expression fl2.24p 


with 


J{iA = E 

n=l 


2n 


e”^Z(n), 


(4.3) 


Z{n) =tre“"-^ 


(4.4) 


is valid only for small e^. To achieve the large p expansion from small e^, in [23] we utilized 
a reciprocal formula which follows from 





£ T cr 


71 

sin Tia 




(4.5) 


This manipulation was further generalized in [33] by using the Mellin-Barnes integration rep¬ 
resentation. Namely, we rewrite the grand potential as an integration 

JM = -l ( 4 . 6 ) 

with 0 < e < 1 and evaluate it in both the regions p > 0 and p < 0. Assuming p < 0, we can 
reproduce the series expansion (14. 3 h by collecting the residues of the integrand in Re(t) > e. 
Assuming /r > 0, on the other hand, we can evaluate the integration by pinching the contour 
so that it encloses the region Re(t) < e. As a result, we obtain the large p expansion of the 
grand potential from the residues of the integrand in Re(t) < e. 
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To explicitly study the grand potential, we use the WKB h-expansion, as in the ABJM 
theory [15]. The h-expansion of Z{n) takes the form 

- OO 

Z(n) = -Y,ti‘‘Zt(n), (4.7) 

^ t =0 

where the overall factor 1/h is due to the normalization by the unit volume of the phase space. 
Note that, since quantum corrections contain h only through ih, the hermiticity of the Hamil¬ 
tonian ensures that quantum corrections only appear in even powers of h. Correspondingly, 
we also decompose the grand potential as 

-| OO 

(4.8) 

^ £=0 

Below we hrst compnte the classical limit Zq^u), by neglecting the ordering of the operators 
and performing the phase space integral explicitly. Then, using the Mellin-Barnes integration 
representation 04.61) we obtain the exact large /i expansion of Jo(h)- After that, we proceed 
to the quantnm ^-corrections and determine J 2 (h) by the same method. 


4.1 Classical limit 


In the classical limit, all the operators can be regarded as c-nnmbers and the trace is the 
(Q, P)-phase space integral divided by Anh. As a result, Zo{n) is 


Mn) = j 

with the classical Hamiltonian Hq given by 

Ho={r- 2)U + (r - 2)T - 2S. 


(4.9) 


(4.10) 


Here f/, T and S are given by 04. 2 p with the operators Q and P replaced simply by c-numbers 
Q and P respectively. Then, the integration in 04.9p is fonnd to factorize as 


Zn(n) = 


{2n)\ 

Att 




' f dQ (sinh^)^“ 


' f dP (sinh^)^^ 

[j (2a)!(2cosh|)’'"-“J 


[J (26)! (2 cosh 


a,6>0,a+6=n 

Using the integration formnla (a G Z>o) 

7“ dx (3mh|)^“ ^ r(^-a)r(f) 
(2a)! (2cosh|)"* 2'^ ■ a]T{m) 


(4.11) 


(4.12) 
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which can be derived recursively by integration by parts, starting with 

1 _ r(f)" 


dx- 


and the formula 


E 


(2 cosh l)™^ r(m) 

r(f-a)rw r(|-t)rM _ r(5±s) 


(4.13) 


(4.14) 


“!r(f)r(x - 2a) «r(f)rt!/ - 2b) 2-2- . - n )' 

which can be shown by considering the generating function with respect to n, we hnally obtain 
the following expression for Zo{n) 

,2 ^\2 


Zo{n) = 


1 r(2n + l)r((|-l)n)'r(|n)' 
dvr r(n + l) r((r — l)n) r(rn) 


(4,15) 


Plugging this into the Mellin-Barnes representation fl4.6p and collecting the residues in 
Re(t) < 0, we obtain the exact large /i expansion of the classical grand potential 




(4.16) 


Here the hrst three perturbative terms come from the residue at t = 0. The coefficients Cq 
and Bq are consistent with the classical Fermi surface analysis in section [3l and the constant 
Ho is 

A„ = h2)(i + JE), (4,17) 

TT \r r — 2 J 

The non-perturbative part consists of three kinds of instantons 


= E 


^ ( 1 ) 

= l^cye r 

e=i 




m=l 


E 

n=l 




(4.18) 


with 


= 


( 1 ) _ (20! r(f)r(-f)r(-<iyH)^ 

7rr(H)2 

1 p / 2m ( 4m ^ "P ^ 

^ \r-2)^ \ r- 2 )^ \ r- 2 ) 


= 


7 r(r - 2)2(m!)2 p^_ 2 (r-pm ^p^_ 2 rT|^ 


2m 


Cm _ I 

^ - 2 


+ 2^/1 


4m 


- (r - l)i}( - 


r — 2 
2{r — l)m 


+ (r — 2)-^(m + 1) + 


rm 


— r'0( — 


2rm 


c(3) = 


r — 2 

(-!)-■ r(|)r(-ii^)¥(-g) 
8?™! r(-^^^)r(-f) 


(4.19) 


where 'ijj{x) is the di-gamma function 'i/j{x) = (92;logr(a:) 
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4.2 Quantum corrections 


Now we shall go on to the quantum corrections. As in [15], with the help of the Wigner 
transformation 


(-^)w 


271 Y 2 


X 



iQ'P 

e 2 « , 


(4.20) 


the trace of operators in Z{n) can be expressed as an integration of a c-function 

Z{n) = j ( 4 . 21 ) 

Practically, the Wigner transformation can be computed by 


f{Q)w = /(Q), f{P)w = f{P), {XYY = {X)w * (?)w, (4.22) 

where the star product is given as = exp[zfi( d q d p — d p'^g)]. In this formulation, we can 
compute all the ^-corrections systematically through the i*r-product. 

As in the case of the A-type theories with J\f = A supersymmetry, there are two sources of 
fl-corrections: the deviation of i^w from Hq, and the deviation of (e“”"^)w from . The 

latter can be incorporated in the same way as in the A-type theories, by applying the Taylor 
expansion 


f{X)w - ^ Qi{XY - ((X - Xw)^)^, (4.23) 

with /(x) = The former deviation can also be calculated similarly as in the A-type theo¬ 

ries, by the Baker-Campbell-Hausdorff formula and the i*r-product. In deriving the Hamiltonian 
operator, all we have to do is to utilize the formula in appendix A of [I5] twice, combining 
first and then to in fl4.ip . Note that S' is a composite of Q and P in this 

case. This again can be treated by the formula fl4.23p with /(x) = logx. In summary, Z{n) 
is given by 


/ 47r^ 




e=i 


-n) 

i\ 


■{Hw-Ho 




e =2 


[-n) 




(4.24) 


Though the calculation is now rather straightforward, let us note that the expression is 
simplified if we introduce 

O P O + P 

t /2 = 2cosh—, T 2 = 2cosh—, S '2 = 2sinh—-—. (4.25) 
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For example, the ^-corrections relevant to Z 2 {n) are given by 




w 


H^ = h^ 


(r — 2)(r^ -|- 2r -|- 2) ^ (r — 2)(r^ -|- 2r -|- 8 ) (r — 2)(r — 4) 


12 t/| 


24T| 


12^1 


(r — 2 )^(r-|-4) (r — 2 )(r-|- 6 ) 2 r(r — 2 ) 


417' 


4T' 


6 t/|T| 




3TiSl 


3 T 2 SI 3U2S^ 


202 


2(^2 _ 3 ^ + 8)U^T^ 2(r - 2 )(r + 2 )t/'^ ^ “ 2 )( 2 r + 


^2(-17w) — 
GsiHw) = 


3 U 2 T 2 SI 

r 


2-L2^2 
2 2 r 


2 r 


17|T| 1/2^2 T2,S2 


2^2 


-2 ^2 


r(r2-|-4) r(r2-|-4) 


4t/2 


4T2 


3U2TiS2 

+ Oih^), 

r 2 2 r^ 


3 UIT 2 S 2 


+ Oih^), 


^ 2 r(r ~ 2 ) ^ 2 r(r — 2 ) 






4r-2f/'T^ 


4r2f/',S.' 

[/aTaFf ' [4^2 ' f/|^ 2 ^ 2 j 


+ 


.2^2 ^ 4r2T^F' 


+ 0{h^), 


(4.26) 


while = 0{h‘^) for i > A. After substituting these into fl4.24p . we can integrate the 

resulting expression by the same technique as in Zo{n). We hnally obtain Z 2 {n) as 


r2(r-2)2n2(l-n)(l + 2n) 

^ ^ 96(1 + (r - l)n)(l + rn) ^ ^ 

By a similar, though more lengthy, calculation, we also obtain Z 4 {n) as 

r^(r — 2)2n^(n — l)(2n -|- 1) 


(4.27) 


ZAin) = 


92160(1 + rn){3 + rn){l + {r — l)n)(2 + (r — l)n)(3 + {r — l)n) 

(8 - 5r + r2)(96 + (-110 + 82r)n + (326 - 58r + 17r2)n2 + (92 + 124r - 5r^)n^ + 14r2n^) 

+ (-432 + 226r)n + (-1616 + 530r)n2 + (-928 + 144r)n3 - 56rn^j Zo{n). (4.28) 


Now let US consider the large /i expansion of the quantum corrections to the grand potential. 
Remarkably, both Z 2 {n) and Z 4 {n) are expressed as Z^iji) times some rational function of 
n. Therefore, J 2 (h) 74 (/r) have the same three species of instantons as 7o(/i), since the 

inhnite sequences of poles of r(f)r(— f)^o(l^) in Re(f) < 0 remain unchanged in the quantum 
correct ions 0 On the other hand, the perturbative part which comes from the residue at n = 0 
is 


72 (h) 

74 (h) 


r(r-l) r(r-l)^ 


487r 


-h 

l)(r 2 


247r 
- 5r + 8 ) 


86 407r 


+ 72 ^(h), 

7r(h). 


(4.29) 


** Though the rational functions have a finite number of poles with n < 0, all of them are canceled by the 
zeroes of [r(rn)r((r — l)n)]“^ in Zo{n) and none produce new instanton effects. 
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Combining the classical value fl4.17p and the quantum corrections fl4.29p for the constant part 
A, we hnd that the result is consistent with fll.lOp . at least up to 0{h^). This is reminiscent 
of a similar relation discovered for the {q,p)k models among the theories of the A-type quiver 
in [23]. Note that the Planck constant seems doubled compared with the {q,p)k models. This 
can be understood by comparing two identical quivers and ^ 3 , where our case with uniform 
Sm = 1 is identified with the (3, 1 ) 2 *: model. See figure HI 



Figure 4: Dynkin diagram of D^. The edges can be determined e.g. from the inner product 
of the canonical basis = (1,—l,0),e^/ = (—1, —1, 0), Ci, = {0,1, = (0,1,1). In 

our parametrization (ll.3p . the levels are given by {kyi,k^,ku,k^i) = 2fc(l, 0, 0, —1). This is 
nothing but the Dynkin diagram for the (3, l) 2 fc model, whose array of Sm is given by 
(v, Sf,, s^, s^/) = (1,1,1, -1). See figure 1 in [23|. 



5 Summary and discussion 

In this paper we have studied the partition function of the superconformal Chern-Simons 
theories of the D-type quiver, and have shown that we can rewrite the partition function into 
that of the Fermi gas system as in the case of the A-type quiver. We find that, again, the 
perturbative corrections of the partition function are summed up to the Airy function, if the 
Hamiltonian of the Fermi gas system is hermitian. Though, for the general D-type quiver, 
in section [3] we only consider the perturbative coefficients in the classical limit fc —)■ 0, the 
Fermi gas formalism is very powerful and allows us in principle to determine the quantum 
corrections and the non-perturbative instanton corrections. 

To further proceed to studying the membrane instanton of the general D-type quivers 
quantum-mechanically by the WKB expansion, it is, however, difficult to handle the non- 
commutative operators in the density matrix, or the exact integration over the phase space 
without taking the large p limit. In the theories of the general A-type quivers, we have 
overcome the difficulties [23II25] by restricting ourselves to those with AA = 4 supersymmetry 
Similarly, here in the theories of the D-type quivers, the difficulty is resolved by choosing 


18 




















the quivers with uniform as in section SJ For these special values of the levels, we have 
found that the non-perturbative corrections consist of three kinds of instantons, and have also 
observed that the constant A can be expressed in terms of that in the ABJM theory (at least 
up to 0(k^)). These are reminiscent of the results for the theories of the A-type quivers with 
the Af = 4 supersymmetry [231121] . 

It is interesting to see whether the symmetry is enhanced for these cases with uniform 
Sm- Also, we hope to interpret these instanton exponents from the dual supergravity picture, 
as membranes wrapping on the tri-Sasaki Einstein manifold, though the geometry is more 
complicated than that for the A-type quivers. Furthermore, we hope to proceed to all the non- 
perturbative corrections including the worldsheet instantons which have not been discussed 
at all in this work. 

After seeing the Fermi gas formalism for the theories of the A-type and D-type quivers, it 
should be interesting to ask whether a Fermi gas formalism exists also for the E-type quivers. 
Also, it is interesting to study other quivers with orthosymplectic groups in [9] from the Fermi 
gas formalism. See e.g. [2B] . 


A A pfaffian formula 


Proposition. Let (0a)i<a<2V and ('^ 6 )i<f,< 2 Ar be functions on a measurable space. Then we 
have 

nv , . 

—^det ( (0a(a:i))l<a<2iV (V'a(a:i))l<a< 2 Jv) = (-l) 2 (^-l)Afpf P (A.l) 

JVl V l<i<N l<i<N J ^ f U.U1 V y 

with the skew-symmetric matrix P 



Pab — 


Dx{(l)a{x)'ipb{x) - 0b(a;)V’a(a;))- 


(A.2) 


Remark. The dehnition of the pfaffian for a skew-symmetric matrix P is given by 

1 ^ 

ptP= ^ (A.3) 

crGS2N ^—1 

Proof. We can prove it by skew-symmetrizing the matrix elements, 

^ I {{4>a{Xi))2NxN {'^a{Xi))2NxN^ 

1 r ^ 

(J&S 2 N *=1 
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(A.4) 



B Another pfaffian formula 


Proposition. Let Pab{x, y) with a, & = 1, 2 be functions of two variables satisfying Ptaiu, x) = 
— Pab{x,y). Let P be a 2N x 2N skew-symmetric matrix consisting of four N x N blocks Pab 
whose (i,j)-component is Pab{xi,Xj). Then, we have 

ptP = \/det(I - zTiP), (B.l) 

N=0 

with various matrices on the right-hand side dehned by the identity operator as 


n = 





(B.2) 


Here the pfaffian on the left-hand side is the finite dimensional one, while on the right-hand side 
the determinant denotes simultaneously the 2x2 determinant and the Fredholm determinant. 

Remark. This is the continuum limit N^o —cxd of the following proposition (See e.g. Propo¬ 
sition 2.1 in [35]). Note that, in taking the limit, we use pf(r2 -|- zP)"^ = det(f2 -|- zP) = 
det(/ — zD.P), which follows from hi ^ = —hi and det hi = 1, and £x the overall signs by 
setting P to be zero. 

Proposition. Let Pa,b with a, & = 1, • • • , 2N^ be a skew-symmetric matrix. Then, we have 

iVoo 

(-1)5(^--I)^“pf(n + ^P) = pfP{4, (B.3) 

N=0 l<si<---<Sf^<Noo 


where Pj^} consists of four N x N blocks whose (i, j)-component is given by Ps^(+N),sj{+N)- 



( {Psi,Sj)NxN 
\{Psi+N,Sj)NxN 


{Psi,Sj+N)NxN 

iPsi+N,Sj+N)NxN 


(B.4) 
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